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Abstract 

We employ genetic algorithms (GA), which allow for an unbiased search for the global minimum 
of energy landscapes, to identify the ordered equilibrium configurations formed by binary dipolar 
systems confined on a plane. A large variety of arrangements is identified, the complexity of which 
grows with increasing asymmetry between the two components and with growing concentration of 
the small particles. The effects of the density are briefly discussed and a comparison with results 
obtained via conventional lattice-sum minimization is presented. Our results can be confirmed by 
experiments involving Langmuir monolayers of polystyrene dipolar spheres or superparamagnetic 
colloids confined on the air-water interface and polarized by an external, perpendicular magnetic 
field. 
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Investigations of the structural and thermodynamic properties of colloids confined at fluid 
interfaces represent a very active topic of current research. Unlike particles in the bulk, the 
effective interactions between the colloidal particles at fluid interfaces are influenced not 
only by the properties of the particles and the solvent but, in addition, by the surface and 
line tensions of the interface.^ Here we focus on the self-assembly scenarios of binary dipolar 
colloids, a system for which various experimental realizations exist, and whose interactions 
have been quantitatively established. Experimentally, studies of two-dimensional ordered 
arrangements of colloidal particles can be fairly easily realized by investigating polystyrene 
particles floating on an oil- water interface^ (system I). Their size, typically lying in the mi- 
crometer domain, allows a direct observation of the particles in light microscopes. Computer 
simulations for the binary casej^ mimicking the experimental setup sketched in Fig. [TJa), 
revealed a surprisingly rich spectrum of exotic ordered equilibrium structures. Another re- 
alization of such systems is offered by the setup of Maret and coworkers, 1 5 1 6 1 7 1 8 1 9 1 employing 
superparamagnetic colloids suspended on a pendant water droplet, as sketched in Fig. [IJb) 
(system II). Here, a magnetic field B applied in the perpendicular direction polarizes the 
particles and induces a repulsive interaction among them that scales with interparticle sep- 
aration r as ~ r~ 3 , see Refs. l9fT0l 

For system I, Sun and Stirner^ have derived the following expressions for the pair po- 
tentials, $jj(r), acting between the two species of particles: 



Here, Ri and Pi are the radius and the dipolar moment of species % and e is the dielectric 
constant of water. In what follows, we assume that Pi = aRf, with some proportionality 
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FIG. 1: (a) Schematic representation of system I (inspired by Ref. H|: two particles of different 
size floating at an oil-water-interface. The dipole moments P\ and P2 are given by the vector sum 
of the dipole moments on the particle-oil-interface, since the dipole interaction is screened in the 
aqueous phase, (b) Schematic representation of system II: superparamagnetic colloidal particles 
trapped on the water/air interface in a pendant water droplet. The external magnetic field B is 
used to tune the interactions between the spheres. 

constant a. Introducing z = Rb/Ra < 1 and x = r/d^ (g?a = 2-Ra), and factoring out a 
prefactor common to all three interactions $jj(r), we arrive at the following expressions for 
the dimensionless interaction potentials ^(x): 
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The functions ^ij(x) are displayed in Fig. |2](a) for two different values of z. Expanding 
the logarithms for x ^> 1 yields the approximate expressions ^ij(x) = (ziZj) 5 ^ 2 /x 3 , where 



1 or z if i,j = A or B. Fig. [27b) shows the asymptotic, power-law form of the 



interactions. In two dimensions, where the exponent of the power-law exceeds the dimension 
of space, the subtle issues of the shape-dependence of the thermodynamics, which arise for 



three-dimensional dipolesp^ are absent. 

For system II, on the other hand, where the superparamagnetic colloids have a suscepti- 
bility ratio m < 1, the dimensionless interaction potentials read 9 * 10 * 1 -^ as ^ij(x) = mi-mj/x 3 , 
where Zij = 1 or m if i, j = A or B. If the overall density of system I is sufficiently low, so 
that the asymptotic forms of ^ij(x) above hold, systems I and II become equivalent, with 
the correspondence z 5//2 m. The ground states of system II have been recently analyzed 
with conventional methods in the work of Assoud et al™ for 0.08 < m < 1, being inspired by 
the phase diagram of a two dimensional, binary hard-sphere mixture.^ Here, we concentrate 
on system I in the domain 0.1 < z < 0.5, corresponding to much stronger asymmetries (i.e., 
3 x 10~ 3 < m < 0.18) but still having sufficient overlap in the parameter range with that 
of Ref . HH to allow for quantitative comparisons. Finally, the system is characterized by the 
number-density rj or its dimensionless counterpart, r]d\. 




FIG. 2: (a) The interaction potentials of Eq. ([TJ for z = 0.3296 (full lines) and z = 0.5 (broken 
lines), (b) Double logarithmic plots of the ^ij(x) for z = 0.5, demonstrating their ~ x~ 3 
power-law dependence for large x- values. 

It is well-known that binary mixtures tend to show a board spectrum of rather complex 
alloy phases, most of them being extremely hard to guess. With the ambition to cover as 
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large a variety of structures as possible and in an effort to cope with this problem in an 
efficient way, we employed genetic algorithms^ (GAs) instead of following the conventional 
approach to the problem of finding minimum energy configurations (MECs), which relies 
on preselected sets of candidate structures. GAs were invented in the 1970s by Holland 
et al. to solve high dimensional and complex problems in engineering science.^ They are 
optimization techniques, modeled after the natural process of evolution, and mimic certain 
biological mechanisms, such as mating and mutation, to find the optimal solution to a 
proposed problem. Due to its special design, a GA is able to take the whole search space 
into account at once and at the same time to concentrate its computing efforts on promising 
regions. It is this global scope that makes GAs an efficient and widespread tool in fields 
like economics and engineering. In the realm of structure optimization, a typical problem 
in atomic condensed matter physics, they have found applications in determining optimal 
atomic clusters,^^ and, very recently, in optimizing extended spatial structures. 1 19 1 20 1 In soft 
matter, applications of GAs to predict equilibrium crystal structures for one-component 
systems has already delivered remarkable results, both in twcP^ and three dimensions] 22 * 23 * 2 ^ 
but binary systems have not been looked upon with GAs so far. 

Based on the approaches presented in Ref. [25j we designed a method to determine ordered 
MECs for binary monolayers. In our method, lattice parameters are freely optimized with 
respect to the free energy, which, at T = 0, reduces to the lattice sum U of the ordered 
structure. The efficiency of the GA allows us to perform our search for equilibrium structures 
among all possible lattices, without posing any bias on the algorithm whatsoever. The only 
limiting factor constraining our search is the maximum of particles per cell that the algorithm 
can handle in reasonable time. The risk of overlooking relevant structures at any point in 
the process is thus minimized. 
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In general, we have to deal with non-Bravais lattices, with s(> 1) particles per unit 
cell. We parametrized these lattices with ha particles in the unit cell belonging to species 
A and the remaining n-Q = s — nj^ particles to the other species, B, coding the position 
vectors of all particles in a binary fashion. We obtained results for particle size ratios in 
the range of z — 0.1 to z — 0.5, including the value of z — 0.3296 that corresponds to 
the particle size ratio used in previous Monte Carlo (MC) simulationsP The concentration 
of small particles, C = n-g/s, was systematically varied for every size ratio. Most of the 
calculations were performed with a maximum of eight particles per unit cell, so C lied in 
the range of 1/8 < C < 7/8. Additional calculations were performed for C = 7/9, as exotic 
phases were to be expected for this particular ratioP^ 

Before presenting the results obtained by the GA, let us discuss some generic expected 
features of the MECs, arising from the functional form of the potential. In sufficiently dilute 
one-component systems, a hexagonal lattice is formed since *&ij(r) ~ r~ 3 at long distances. 
This "default structure" is also found in two limiting cases of the binary mixture: first, for 
z — > 0, when the smaller B-particles are of vanishing size compared to the larger A-particles, 
and second for z — > 1, when the two species become indistinguishable. In the first case 
we expect the large particles to form hexagonal lattices, more or less unperturbed by the 
presence of the small ones. If z is small, the B-particles can form independent patterns in the 
interstitial regions of the hexagonal lattice formed by the A-particles without paying a high 
penalty in energy. Thus, distinct groups of small particles should be observed with their 
size fulfilling the stoichiometric requirements. In the opposite the similarity in size 

increases, the B-particles require more space. Now the concentration and the dependence 
of the interspecies interactions on the distance are the key quantities that determine the 
ordered equilibrium structure. 
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FIG. 3: Minimum energy configurations for C < 1/2, i.e. < «a f° r = Particle 
diameters are not drawn to scale. 



In Figs. [3] and [4] we present the ordered equilibrium configurations found with our GA- 
method for increasing concentration of small particles C at an overall particle density of 
r]d\ = 0.6. To facilitate discussion we divide our results into two blocks: structures with 
n B < (Fig- 13| and structures n B > (Fig. [1]). We begin with the case n B < n&. For 
sufficiently large values of z, i.e., z = 0.3296 and z = 0.5, the small particles are able to 
perturb the (ideal) hexagonal structure formed by the large particles in their immediate 
surrounding, even though they represent the minority component. In regions where there 
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FIG. 4: Minimum energy configurations for C > 1/2, i.e. n-Q > for rjd\ = 0.6. Particle 
diameters are not drawn to scale. 



are no small particles in the immediate neighborhood, the hexagonal structure of the large 
particles prevails, whereas the small particles are found in the center of squares formed 
by the large particles [see Figs. |3[i),(j),(m) and (n)]. On the other hand, if z is below 
a certain threshold, i.e., z = 0.1 and z = 0.2, the influence of the small particles is not 



sufficient to cause a substantial modification of the hexagonal pattern of the large particles. 
Instead, since the small particles experience mutually a very weak repulsion, they tend to 
stay close together in this ^-regime, arranging in lanes which meander through the hexagonal 
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lattice formed by the large particles, see Figs. |3^c), [3^d), |3](g), and[3](h). For very small 
concentrations, i.e., C = 1/7 and C = 1/5 [cf. Figs. [3^a), |3](b), |3^e), and[3^f)], the 'lanes' of 
the small particles are interrupted by intervening big ones due to the simple fact that not 
sufficiently many B-particles are available in the system. It might be possible that also in 
these cases pure B-particle lanes will form but this would require a much larger unit cell, 
which was not included in our study. 

When the small particles become the majority component, the found structures become 
much more complex, as can be seen from Fig. |4} As expected, B-particles are observed to 
arrange in distinct groups for many parameter settings. We find small particles forming 
dimers, [Figs. Qb), |3|c), and[4|d)], elongated- [Fig. [l|g)] and triangular-trimers [Figs. |4](c) 
and pjje)], as well as chain-like pentamers [Fig. ml)] and heptamers [Fig. [4^h.)] . In con- 
trast, the A-particles form rather simple lattices which accommodate in their interstitial 
regions these sometimes rather complex groups of B-particles. For C = 2/3, dimers of small 
particles, observed for moderate values of z [Fig. [4](k)] are a precursor of lane formation 
[Fig. |4^p)]. For small values of z, the large particles form a hexagonal pattern and the small 
particles are distributed in the interstitials. At z — 0.3296 the hexagonal structure of the 
large particles is distorted and the small particles are grouped in dimers. If the particle 
size ratio is further increased, the dimers of small particles change their orientation and 
lane formation sets in which now strongly distorts the underlying lattice of A-particles. The 
scenario repeats itself for C = 5/7, [vertically from Fig. mb) to Fig. |4](q)] but here the stoi- 
chiometry, which does not accommodate a single B-particle in the interstitials of A, forces 
a much richer structure: for z = 0.1, Fig. |4](b), the B-particles form an ordered array of 
monomers and dimers, which transforms into an array of monomers and linear trimers for 
z = 0.2, Fig. |4^g). For z = 0.3296, the aggregates of B-particles become monodisperse, 
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zig-zag-like pentamers, Fig. [4^1); and, finally, for z = 0.5, formation of two interchanging 
kinds of B-lanes, thick and thin, takes place, Fig. |4][q_) . 

Evidently, the higher the value of C, the more surprises are hidden in the system. Let 
us discuss the sequences for C = 7/9, C = 4/5, and C = 6/7 in more detail. For C = 6/7 
and for small ^-values, the interstitials of the A-lattice are occupied by three different kinds 
of B-aggregates: monomers, dimers and triangular trimers, Fig. ^c). For z = 0.2, B- 
particles arrange in zigzag-shaped heptamers which distort the hexagonal pattern of the 
large particles, Fig. |4]^h). However, in contrast to the smaller C-values a further increase 
of z does not directly lead to lane formation but rather a new, exotic structure intervenes 
for z = 0.3296, Fig. |4^m). Here, two neighboring heptamers merge, forming thereby a 
sequence of alternately oriented cup-like structures, each of them hosting an A-particle. 
If z is increased further, thick and thin lanes are again formed, Fig. |4jj). In the second 
sequence, C = 4/5, the B-dimers [Fig. ptd)], observed for z = 0.1, transform into zig- 
zag lanes [Fig. |4^i)], pearl-necklace- lanes [Fig. |4](n)] and finally into rings, each of them 
surrounding one A-particle, [Fig. |4Ts)], as z grows. Finally, for the sequence C = 6/7, a 
ring-like structure is formed for z = 0.1 to z — 0.3296, Figs. [4|V), |2]Q) and|4^o), where every 
large particle is surrounded by six triangular trimers of small particles, forming a structure 
resembling a Kagome lattice [Fig. Qe)]. The B-interparticle distance within the trimers 
increases with z. Finally, for z = 0.5 lane formation sets in once more [Fig. |4](t)], but now 
the lanes formed by the small particles are interconnected, due to the high concentration of 
B-particles. 

The structures found by the GA for z = 0.5 correspond to m — 0.177 in the terminology 
of Assoud et al. 14 For those stoichiometries that have been considered both in the present 
work and in Ref. HH the same structures were found by both the conventional approach and 
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FIG. 5: Two sequences of MECs for z = 0.5 and two different concentrations C. Structural change 
is observed in both sequences, once between rjd\ = 0.6 and rjd\ = 0.8 (top) and once between 
■qd\ = 0.4 and r/d\ = 0.6 (bottom). 

by the GAs. The GAs offer a higher predictive power and flexibility to identify structures in 
the regime of large size-disparity and high C-values, where the unit cells become increasingly 
complex and conventional methods reach their limitations] 14 * 26 ^ 

Were the interaction potentials of Eq. ([TJ to be pure power-laws, as those employed in 
Ref. H3J the overall density r]d\ would be an irrelevant parameter, in view of the absence of 
any length scale in the interactions. Since the logarithmic dependence, as well as the hard 
cores of the present interactions, Eq. ([lp, set in at very small separations [cf. Fig. |2|b)], 
we expect that the structures reported in Figs. [3] and [4] will be stable for a broad range of 
small to intermediate densities. To put this assumption into a check, we have investigated 
the ordered equilibrium configurations of our system varying the density rjd\ from 0.2 to 
0.8. We report selected results in Fig. [5] For concentration C = 1/3, upper row, the lane- 
like structure remains unchanged up to a density at least as high as r)d\ = 0.6, only the 
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distances are scaled as r\ increases. As the density is further increased to r\d\ = 0.8, the 
following structural transition takes place: every second B-lane is dissolved by merging one 
of the neighboring ones; consequently, the inter-lane distance and the density in each of these 
lanes is increased by a factor of two. For C = 5/6, a structural change induced by density 
sets in earlier, namely between rjd\ = 0.4 and r\d\ = 0.6. Here, the original pattern shows 
rings of B-particles while, again, at high densities lane formation is observed. Moreover, all 
structures shown here have the lowest energy values found by the GA but, at the same time, 
other ones with very small energy differences have been discovered in some GA runs. Both 
the effects of the density and the discussion of 'quasi-degenerate' ground states will be the 
subject of a future publication. 

We have applied genetic algorithms to examine ordered equilibrium configurations for 
binary, two-dimensional dipolar mixtures. Despite the simplicity of the interactions involved, 
the system shows a tremendous variety and complexity of the resulting structures. This 
convincingly demonstrates the power of this novel optimization technique to deal with the 
problem of finding the global minimum of a rugged potential energy surface, which becomes 
increasingly involved as the number of components of the mixture increases. We have 
demonstrated that GAs are an effective and reliable new tool which merits a more widespread 
appreciation in the soft matter community, on equal footing with complementary tools, such 
as Monte Carlo simulated annealing and Molecular Dynamics simulations. The variety of 
the encountered structures can easily be verified experimentally considering either a binary 
mixture of two-dimensional superparamagnetic colloidal particles^ 21 or a binary mixture of 
polystyrene spheres floating on an oil-water interface^. 
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